September 27, 2012 



Lattice results on charmonium 



Daniel Mohler 
TRIUMF, 4 004 Wesbrook Mall Vancouver, BC V6T 2A3, Canada 



PRESENTED AT 



The 5th International Workshop on Charm Physics 

(Charm 2012). 



1 Introduction 



Lattice QCD (LQCD) is a tool for ah initio calculations of the theory of quantum 
chromodynamics (QCD) in the non-perturbative regime. As such it is ideally suited 
for calculations of the hadron spectrum and to calculate properties of hadrons. In the 
charmonium spectrum, many observed states above multiparticle threshold do not fit 
into the picture resulting from simple potential models such as [TJ. While these states 
can be studied with a variety of models, QCD calculations are very desirable. Fur- 
thermore, the rich spectrum observed in charmonium is an ideal benchmark for lattice 
calculations. The strange and light quark masses and the lattice scale are usually de- 
termined from light-quark observables. This makes calculations of the charmonium 
spectrum an ideal benchmark calculation. 

While calculations in LQCD were limited to larger than physical pion masses for 
many years, this obstacle has been overcome in recent years and there have been first 
calculations at (2j |3] or close to jl] physical pion masses using at least 2 flavors of 
dynamical quarks. 

As an example, some ground state masses of light-quark mesons and baryons have 
been calculated with full control of systematic uncertainties [5]. In these proceedings I 
will highlight corresponding progress for the calculation of the charmonium spectrum 
from Lattice QCD. In Section [5] some basic issues for lattice calculations are briefly 
reviewed. Apart from discussing continuum, infinite volume and chiral extrapolations 
a well-established method to extract excited states on the lattice is presented. In 
section [3] results from various recent simulations are collected. Section H] provides 
some brief concluding remarks. 

2 Lattice methodology 

2.1 A threefold of extrapolations 

In LQCD observables are calculated using Monte-Carlo methods. Any such simulation 
on a finite number of configurations, in our case the finite number of dynamic gauge 
configurations within an ensemble, will have a stochastic uncertainty associated with 
the estimate. In addition to the stochastic uncertainty, rigorous quantitative com- 
parisons with experimental results can only be made after systematic extrapolation 
to the continuum and infinite volume limits and after extrapolation or interpolation 
to physical quark masses. How important these uncertainties are depends on the 
observable and on the lattice action used. In this section, we discuss the role of 
these extrapolations and provide a short explanation how these issues are handled in 
practice. 

While a fully systematic treatment and a quantitative comparison to experiment 
mandate these extrapolations, many qualitative insights can be obtained from sim- 



1 



ulations where some of the uncertainties associated with continuum, infinite volume 
or chiral extrapolations have not been fully controlled. However, caution has to be 
exercised when interpreting the resulting data. In Section |3] we will highlight this us- 
ing the example of the charmonium hyperfine splitting, which can be quite sensitive 
to discretization effects. 

2.1.1 Continuum limit 

Observables measured on the lattice contain discretization errors which depend on 
the lattice spacing a. To compare to experiment an extrapolation to the continuum 
has to be performed, which corresponds to the limit a({3,m qi ) — > 0, with the inverse 
gauge coupling j3 — ^ and the quark masses m qi . Given multiple lattices of the same 
physical size with different lattice spacings one can fit the data with the functional 
form expected for the leading discretization errors and obtain values in the continuum 
limit. In practice this is a delicate business in simulations with dynamical fermions. 
In particular the extrapolation should be performed along lines of constant physics 
and the physical volume is usually not constant for the initial estimate of simulation 
parameters. It should be stressed, that while complicated in detail, this is a straight- 
forward procedure for many observables. 

2.1.2 Infinite volume limit 

For a given volume, interactions due to the periodicity of the lattice lead to finite 
volume effects. Provided a large enough lattice, these cause exponentially small cor- 
rections to the spectrum of hadrons [36] . These are usually given by the lightest 
particle in the spectrum, the pion, and are therefore of order exp(— m n L), where m n 
is the pion mass and L is the size of the box. It is often assumed that these effects 
can be neglected for lattices where the dimensionless quantity m^L > 4. A better 
approach consists of measuring all observables on boxes of increasing size (keeping all 
other simulation parameters constant) and performing an extrapolation to the infinite 
volume limit. 

2.1.3 Physical quark masses 

While calculation at physical quark masses are now possible, present simulations still 
commonly employ larger than physical light-quark masses. In this case an extrap- 
olation to physical quark masses has to be performed. This is often accomplished 
by taking a fit form motivated by chiral effective field theories (for example Chiral 
Perturbation theory (%PT) in the case of light mesons or Heavy Meson Chiral Per- 
turbation Theory in the case of heavy-light mesons). In the lattice literature this 
extrapolation is often referred to as "chiral extrapolation", even in cases where the 
actual extrapolation is to non- vanishing light-quark masses. 
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2.2 Excited states and the variational method 



On the lattice correlators in Euclidean space-time are measured. It can be shown [7] 
that these behave as 

(o 2 (t)Oi(0)^ oc J2 e ~ tE " < °\02\n >< n|Oi|0 > . (1) 

The left hand side can be expressed as a path integral and its discretization can 
be used to calculate this expression using Monte-Carlo techniques. In the simplest 
example, interpolating field operators 0\ and O2 that create and destroy a hadron 
with specific flavor and quantum numbers can be used to obtain information about 
the hadron spectrum. As an example, the interpolator d^u could be used for a ir + . 

Using interpolators of definite quantum numbers it is therefore straight-forward to 
calculate the ground state in each channel, which can be obtained from the large time 
behavior of the Euclidean correlator. To extract not only the ground states but the 
low-lying spectrum one commonly calculates a matrix of correlators at some source 
time slice t{ and for every sink time slice tf 

Cij (t =t f - U) = J2(m(t f )O](U)\0) (2) 

u 

= £ e -^<0|O^>(n|Oj|0> . 

n 

For charmonium interpolators of definite J PC , where J is spin, P is parity and C is 
charge conjugation would be usecfl The low-lying spectrum is extracted by solving 
the generalized eigenvalue problem 

C(t)fW = \ {k \t)C{t )^ , (3) 
A( fc )(t)oce"^ (l + (D(e- tAE *)) . 

for each time slice. At large time separation only a single state contributes to each 
eigenvalue. This procedure is known as the variational method [SI [91 [10]. It is used 
for several of the results presented in Section [3j 



*In practice this is complicated by the fact that, due to the loss of rotational symmetry, an infinite 
tower of continuum states contributes to a given irreducible representation. 
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3 Recent charmonium simulations 

3.1 Results from FNAL/MILC 



Study 


sea 


Charm 


Ops 


States 


Comment 


2009 FNAL/MILC 


2 + 1 


FNAL 


2 


low S, P 


a > 0.09 fm 


2012 FNAL/MILC 


2 + 1 


FNAL 


many 


many 


a > 0.045 fm 


Current FNAL/MILC 


2 + 1 + 1 


HISQ 


2 


low S, P 


a > 0.06 fm 


Future FNAL/MILC 


2 + 1 + 1 


HISQ 


many 


many 


a > 0.06 fm 



Table 1: Analysis campaigns within the Fermilab Lattice and MILC collaboration 
charmonium project. 2009 results refer to [IT] . In addition to the analysis campaigns 
using lattices with 2+1 dynamical flavors of asqtad quarks, there are also current and 
future campaigns using gauge configurations with 2+1+1 flavors of highly improved 
staggered quarks (HISQ), where the charm quark is also included in the sea. 
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Figure 1: Preliminary results for the hyperfine splitting Mj/^, — M Vc as a function of 
the lattice spacing squared. For further details please refer to the text. 

In this section some preliminary results obtained by the Fermilab Lattice and 
MILC collaborations are discussed. The goal of their charmonium project is a pre- 
cision study of the low-lying charmonium states also including higher spin states 
including a full error budget for uncertainties stemming from the Monte Carlo analy- 
sis and for systematic uncertainties. Tabled] provides an overview of past and current 
analysis campaigns related to charmonium. Results presented in this publication are 
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preliminary results from the 2012 analysis campaign using 2+1 flavor asqtad config- 
urations generated by the MILC collaboration [121 E3] • 

Figure [T] shows the lS-hyperfine splitting Mj/^ — M Vc as a function of the lat- 
tice spacing squared compared to experiment. Five different lattice spacings ranging 
from a rather coarse a ~ 0.15fm to a very fine a ~ 0.045fm are used for two differ- 
ent masses of light sea quarks mi, specified as a fraction of the strange quark mass 
m s . Note that both light quark masses are heavier than physical, necessitating an 
extrapolation. Experiment values from the PDG [H] and a recent result by the BES 
collaboration [15] are indicated by magenta and blue stars respectively. While the 
plot does not yet show a continuum or chiral extrapolation, which will be performed 
on the final dataset, the data is consistent with recent experimental determinations. 
Within the renormalization scheme used, there is a small but significant sea quark 
mass dependence of the results. From Figured] it is also evident, that the use of an 
improved heavy-quark action, the so-called Fermilab action [TH] leads to rather mild 
discretization effects. While the current results are fully consistent with experiment 
we would like to point out that contributions from quark annihilation have not yet 
been included and that the physical t] c has a non-negligible width, which is neglected 
in this simulation. The effect of the quark-annihilation contributions will be discussed 
below. 
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Figure 2: Preliminary results for the IS - IP splitting as a function of the lattice 
spacing squared. For further details please refer to the text 

In addition to spin-dependent quantities like the hyperfine splitting, one can also 
take a look at spin-independent splittings. As an example one can define mass cen- 
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troids for the IS* and IP states 

M T s = ^(M Vc + 3M J/ v), (4) 
M TF =^(M Xc0 +3M Xcl +5M Xc2 ). 

Figure [2] shows the resulting 15 - IP splitting as a function of the lattice spacing 
squared. The lattice data for two different sea quark masses has been plotted as black 
and red diamonds and the experimental data from the PDG is indicated by a magenta 
star. Again, within the renormalization scheme used the data show a non- negligible 
sea-quark mass dependence at fine lattice spacings. While extrapolations have not 
yet been performed the results are consistent with the experimental value. 



3.2 Charmonium hyperfine splitting from HPQCD 




Figure 3: Results for the hyperfine splitting from HISQ charm quarks by the HPQCD 
collaboration. Plot taken from [17]. For further details please refer to the text. 

Figure [3] shows results from the HPQCD-collaboration [17] for the charmonium 
hyperfine splittingQ. For the simulation highly improved staggered quarks (HISQ) are 
used for the valence charm quarks and the sea contains 2+1 light flavors of asqtad 
staggered quarks. The simulation data is displayed as blue circles, along with a 
continuum extrapolation of the results (red line). The final results are indicated by 

tin the meantime HPQCD has updated their calculation and the results can now be found in 

[351. 
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the red band and the numerical value obtained is 118.8 ± 2.4MeV. For comparison 
the PDG value for the hyperfme splitting is displayed as a black circle. 



3.3 Disconnected contributions to the hyperfme splitting 

Before presenting further results, it is worthwhile to briefly discuss the role of dis- 
connected contributions to the hyperfme splitting. Their determination is extremely 
challenging and therefore these contributions are often neglected or estimated from 
theory. For qq mesons the relevant quark-line diagrams are displayed in Figure HI Ne- 
glecting annihilation effects corresponds to omitting the quark line diagram displayed 
on the right side of the figure. 



Figure 4: Illustration of quark-line diagrams contributing to Charmonium correla- 
tors. The calculation of contributions from the diagram on the right is technically 
demanding and these contributions should be small. Therefore this contribution is 
often omitted. 

For the results of the hyperfme splitting from HPQCD presented in the previous 
subsection, contributions from disconnected diagrams are estimated from perturba- 
tion theory which leads to a shift of simulation results enlarging the hyperfme split- 
ting. Correspondingly the central value red band in Figure [3] with the final results 
indicates a hyperfine-splitting larger than the continuum extrapolated value. Another 
strategy consist of estimating the effects of annihilation diagrams by direct calcula- 
tion. This approach has been taken by the MILC collaboration in [T8]. The authors 
find that disconnected contributions to the hyperfme splitting reduce the hyperfme 
splitting by l-4MeV. It should be emphasized that novel methods [191 ISHj should lead 
to a more precise determination of annihilation effects in future direct simulations. 

3.4 Two exploratory studies 

In this section we take a look at results from two recent exploratory studies [2Tj 122] . 
While their focus is on heavy-light mesons, the low-lying charmonium spectrum has 
been calculated as well. In contrast to the results discussed in the previous subsection 
the aim of these studies is to push the boundaries with regard to extracting excited 
states and to move towards calculations of hadron resonances on the lattice. As 
these are more difficult problems, current simulations commonly do not include a full 
treatment of systematic uncertainties. In particular they often lack continuum and 
infinite volume extrapolations. 




7 



Nl x N T 


a[fm] 


L[fm] 


#configs 


m^MeV] 


32 3 x 32 


0.0907(13) 


2.9 


> 198 


156 <m n < 702 



Table 2: Parameters of the PACS-CS lattices [I] used in J2JJ. iV £ and N T denote the 
number of lattice sites in spatial and time directions, a denotes the lattice spacing 
and L denotes the physical volume. For the results presented in these proceedings 
the values from the ensemble with the lightest pion mass Mn ~ 156MeV are used. 
Please refer to [21] for further details about the simulation. 



The first study [21] uses Nf = 2 + 1 Wilson-Clover configurations generated by 
the PACS-CS collaboration [I] . For charm , the Fermilab prescription [161 GSj is used. 
In the following, results from the lightest pion mass, which is close to the physical 
pion mass will be presented. For details about the gauge configurations please refer 
to Table H 



Nf x N T 


a[fm] L[fm] 


#configs 


m^MeV] 


16 3 x 32 


0.1239(13) 1.98 


280 


266(3)(3) 



Table 3: Parameters for the Clover- Wilson lattices used in [22]. For an explanation 
of the parameters please refer to Table [21 

The second study uses configurations with Nf = 2 flavors of nHYP-smeared [21] 
Wilson-Clover quarks generated for reweighting studies [251 [26]. F° r the charm quarks 
the Fermilab method [TBI [23] is employed. Details about the gauge configurations are 
listed in Table [3] Further information about this simulation can be found in [22] . 

Figure [5] shows the difference between several low-lying charmonium masses and 
the mass of the spin-averaged IS ground state Mjg = \{M Vc + 3Mj/^). The quantum 
number channels are denoted by either their particle names (for the low-lying states) 
or by their J PC , where J is the Spin, P stands for the parity quantum number and 
C is charge conjugation. For both studies, there is good qualitative agreement with 
the low-lying spectrum observed in experiment and the IS, IP and 2S states are well 
reproduced. In addition to the low-lying states extracted in both simulations, the full 
sets of ID and 2P states are extracted in the more recent study [22]. The high sta- 
tistical precision needed for this has been achieved through the distillation technique 
proposed in [19J. To illustrate the good statistical precision, Table H] summarizes the 
charmonium results. For a full discussion of these results please refer to [22J. 
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Figure 5: Charmonium spectrum as extracted from the lattice from the two different 
studies described in detail in the text. The red crosses (displaces slightly to the left) 
are from [21] while the blue crosses (displaced slightly to the right) display results 
from [22] • The error bars include statistic errors and the systematic uncertainty from 
the scale-setting procedure. Experimentally observed states are plotted as black bars 
or filled boxes. The level corresponding to the well-established X(3872) has been 
plotted for both choices of possible quantum numbers (1 + + or 2 h ). 



3.5 Excited states from anisotropic Wilson-Clover lattices 

The Hadron Spectrum Collaboration has recently studied [20] ground and excited 
charmonium states using anisotropic Wilson-Clover lattices with volumes 16 3 x 128 
and 24 3 x 128. The spatial lattice spacing is a s ~ 0.12fm and the anisotropy between 
the lattice spacing in spatial and time directions is £ = 3.5. The sea-quark pion 
mas for their simulations is = 396MeV. They use the variational method with a 
basis containing interpolating fields with up to three covariant derivatives. A similar 
basis has previously been used for studies of light-quark mesons [2S] • Before turning to 
results let us discuss two notable features about this simulation. First, the distillation 
technique [19] allows calculations with small statistical uncertainty, even for some of 
the higher excitations. In addition advanced spin identification techniques [29] allows 
for unambiguous assignment of certain lattice states to a particular continuum spin. 
In particular, this enables the authors of [2D] to clearly identify Spin 4 states and 
states with spin-exotic quantum numbers. 

Figure |6] shows the results for the charmonium spectrum from [20]. A large number 
of conventional and non-conventional states has been extracted. The states displayed 
as red boxes on the right side of the figure correspond to spin-exotic states which 
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Mass difference 


This study [MeV] 


Experiment [MeV] 


IP — IS 


441.7 ±4.0 ±4.6 


457.5 ±0.3 


2S - IS 


592.3 ±4.9 ±6.2 


606.1 ± 1.0 


IS hyperfine 


107.9 ± 0.3 ± 1.1±q 2 


116.6 ± 1.2 


IP spin-orbit 


39.7 ±2.1 ±0.4 


46.6 ±0.1 


IP tensor 


11.02 ±0.87 ±0.12 


16.25 ±0.07 


IP hyperfine 


3.7 ±2.7 


-0.10 ±0.22 


2S hyperfine 


57.9 ±2.0 


49 ±4 


2P spin-orbit 


24.6 ± 15.7 ±0.3 




2P tensor 


2.2 ±4.3 




2P-1S 


836.4 ±30.5 ±8.8 




271% - il% 


1065.6 


1084.8 ±0.8 



Table 4: Mass differences in the charmonium spectrum in MeV compared to exper- 
imental values (calculated from [T3]; the value for the IP hyperfine splitting is from 
[27]). Bars denote spin-averaged values. For the results of this paper, the first error 
denotes the statistical uncertainty and the second error denotes the uncertainty from 
setting the lattice scale. It is stressed that the gauge ensembles used do not allow 
for a continuum and infinite volume extrapolation. Consequently qualitative but not 
quantitative agreement is expected. In the last line we also provide the splitting 
2Mjj — M-= which can be directly compared to the results quoted by the Fermilab 
lattice and MILC collaborations [IT] and also to the value of 2M^ s — M Vc quoted by 
HPQCD in [28]. 

can be cleanly identified by means of spin-identification [29]. The authors point out 
that some of the non spin-exotic states fit into a possible multiplet of hybrid states. 
For higher-lying states the authors are careful to point out that these states will be 
influenced by unphysical thresholds and should only coincide with physical states 
up to the hadronic width of these states. Notice that the distillation technique [T9] 
allows for the inclusion of these thresholds in a straight-forward manner. For more 
information please refer to [20] . 

To assess the discretization effects present in this study it is instructive to take 
a look at the IS" hyperfine splitting M^fs = 80 ± 1 MeV. Unlike the results from 
the MILC and HPQCD collaborations this result differs significantly from the exper- 
imental value. The difference illustrates the typical size of discretization effects. One 
has to stress that this can be improved by the use of an improved heavy quark action 
or by a continuum extrapolation. We will return to this in subsection 13.61 where it is 
shown that these issues can be overcome by a proper continuum extrapolation. 
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Figure 6: Difference between the mass of charmonium states and the r] c mass M Vc . A 
large number of states including spin-exotic states has been identified. Figure taken 
from [20] 

3.6 Results from the BMW and QCDSF collaborations 



iVf x N T 


a[fm] 


L[fm] 


m^[MeV] 


24 3 x 48/32 3 x 64 
24 3 x 48/32 3 x 64 


0.0795(3) 
0.0795(3) 


1.9/ 2.5 
1.9/ 2.5 


442 
348 



Table 5: Parameters of the QCDSF lattices used in [30] 



In a joint project, the Budapest-Marseille- Wuppertal(BMW) and QCDSF collab- 
orations are investigating the charmonium spectrum using configurations with 2 + 1 
flavors of sea quarks. To this end they use both Wilson-Clover gauge configurations 
generated by BMW [31] and SLiNC gauge configurations [32] generated by QCDSF. 
The BMW configurations span a wide range of lattice spacings 0.054fm < a < 0.092fm 
and pion masses 120MeV < M n < 520MeV on volumes of size L 3 x T = 32 3 x 64 and 
64 3 x 144. For the QCDSF configurations the parameters of used configurations are 
listed in Table |5j For both actions the respective fermion action is also used for the 
charm quarks. 

Figure [7] shows preliminary results for the charmonium spectrum calculated on 
QCDSF configurations compared to the experimental spectrum (with a particular 
assignment for some controversial states). A large number of states including several 
excited states were extracted and the results agree qualitatively with experiment for 
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Figure 7: Preliminary results for the Charmonium spectrum calculated on SLiNC 
configurations [32] generated by the QCDSF collaboration. Figure taken from [30] 

states below the DD threshold. 

Figure [8] shows preliminary results for the charmonium hyperfine splitting cal- 
culated on the HEX-smeared Wilson-Clover gauge configurations generated by the 
BMW collaboration. The results show a strong dependence on the lattice spacing for 
this action. This is in stark contrast to the results employing improved heavy-quark 
actions. Nevertheless a value compatible with the experimental hyperfine splitting 
can be obtained when performing the proper continuum extrapolation 

4 Concluding remarks 

In recent years there has been renewed interest in lattice calculations of the char- 
monium and heavy-light meson spectrum. The availability of configurations with 
physical or close-to-physical sea quarks, improved heavy quark actions and new cal- 
culational techniques like distillation promise many improvements to existing results 
in the near future. While the era of precision calculations using Lattice QCD has 
barely begun, there already is remarkable success. In particular, calculations at or 
around the physical pion mass are becoming common place. There are multiple col- 
laborations generating gauge configurations with 2+1 or 2+1 + 1 dynamical flavors 
suitable for charmonium physics. This rapid progress is also highlighted by several 
more recent charmonium results [331 ES [35] which have appeared since the talk cov- 
ered in these proceedings was delivered. 

Furthermore for states below inelastic thresholds, the extraction of resonance prop- 
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Figure 8: Preliminary results for the charmonium hyperfine splitting calculated on 
the HEX-smeared Wilson-Clover gauge configurations generated by the BMW col- 
laboration. Figure taken from [30] 

erties [361 EH EH [39] from the lattice is within reach. In the charmonium spectrum 
this necessitates the inclusion of DD, DD* and other states and a first step in this 
direction has already been taken [ID]. Nevertheless there is still much work ahead 
with regard to QCD calculations of hadronic excitations. 
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